Susceptibility of the 2D S — \ Heisenberg antiferromagnet with an impurity 
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We use a quantum Monte Carlo method (stochastic series expansion) to study the effects of a 
magnetic or nonmagnetic impurity on the magnetic susceptibility of the two-dimensional Heisen- 
berg antiferromagnet. At low temperatures, we find a log-divergent contribution to the transverse 
susceptibility. We also introduce an effective few-spin model that can quantitatively capture the 
differences between magnetic and nonmagnetic impurities at high and intermediate temperatures. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Cx, 75.40.Mg 



Static impurities can be introduced into the CuC>2 
planes of the high-T c cuprates as a means of probing their 
electronic correlations and excitations; some unwanted 
impurities and defects are also always present in "pure" 
systems. Impurities coupled to a two-dimensional (2D) 
host system hence constitute an important class of quan- 
tum many-body problems. When Cu ions in hole-doped 
CuC>2 planes are substituted with Zn, experiments indi- 
cate that local magnetic moments are induced at the Cu 
sites neighboring the nonmagnetic Zn impurities Q . This 
behavior can be reproduced in spin-gapped Heisenberg 
antiferromagnets |2j, e.g., in systems of weakly coupled 
spin-i ladders [U, y[ Although spin models cannot 
address questions specifically related to a metallic or su- 
perconducting host, static impurities in various gapped 
(paramagnetic) and gapless (antiferromagnetic) Heisen- 
berg systems are important limiting models for under- 
standing the physics of quantum spin defects. They 
also have direct experimental realizations in the cuprates 



Several studies have addressed static vacancies 0, El 
SHU and added spins [H E3 in 2D S = 1/2 Heisen- 
berg models. The localized moments forming in spin- 
gapped systems are now well understood. Recently, a 
universal behavior was predicted for the T > im- 
purity effects in systems that have long-range order at 
T = (i.e., in the "renormalized classical" regime H3)j 
as well as in systems close to a quantum-critical point. 
Some predictions for quantum-critical systems have been 
confirmed numerically [13L whereas other aspects of the 
theory remain unsettled Hfil | . The predictions in the 
renormalized-classical regime have not yet been tested. 

In this Letter, we consider the 2D square-lattice spin- 
i Heisenberg antiferromagnet with (i) a vacancy (non- 
magnetic impurity) or (it) an off-plane added spin (mag- 
netic impurity) coupled to a single host spin. We perform 
quantum Monte Carlo calculations, using the stochastic 
series expansion (SSE) technique jlfij . and determine the 
effects of the two different impurities on the uniform mag- 
netic susceptibility. Our results confirm quantitatively 
a classical-like longitudinal Curie contribution resulting 
from alignment of the impurity moment with the local 
Neel order 0. However, we also find evidence of a loga- 



rithmic divergence as T — > in the transverse component 
of the impurity susceptibility, instead of the predicted T 
independence at low T. We point out qualitative differ- 
ences between the vacancy and the added spin at high 
and intermediate T and explain these in terms of an ef- 
fective few- spin model. 

We begin by defining three spin-i Hamiltonians de- 
scribing the models investigated: 



H a = 
H_ = 

H, = 



J ^ Sj • Sj, 
J ■ Sj, 



<*,3> 



(la) 
(lb) 

(1c) 



where denotes a pair of nearest-neighbor sites on 
a periodic L x L lattice. Ho is the standard Heisenberg 
model. In H- the spin So has been removed, creating 
a vacancy. We also study systems with two vacancies at 
maximum separation on different sublattices. In 7J+ an 
added spin-|, S a , is coupled to one of the spins, So, in 
the plane. We here only consider J± = J. 
We have calculated the total susceptibilities 



x2 = i((5>, 




(2) 



with k = 0, — , + corresponding to Ho, H-, and H + . Wc 
have used square lattices with L = 4, 8, 16, 32, and 64, at 
temperatures down to T = J/32. The number of spins 
N = L? for k — and L 2 ± 1 for k = ±. In order to 
determine the effects of the impurities, we follow Ref. 
and define the impurity susceptibilities 



2,± 

Ximp 



X± - Xo- 



(3) 



In the case of two vacancies, the impurity susceptibility 
is further normalized by a factor i, so that for L — > oo 
it should be the same as for a single vacancy. 

The temperature dependence of xfmp was discussed by 
Sachdev et al. §| on the basis of quantum field theory. 



2 



They concluded that at T = 0, the component Xi mp par- 
allel to the direction of the Neel order vanishes, while 
the perpendicular component Xim P — C3 / ps, where C3 is 
a constant and p s is the spin stiffness of the host. For 
T > the correlation length £ grows exponentially as 
T — > |l2| . For a general impurity spin, a moment of 
magnitude S was then predicted to result from the align- 
ment of S with a large Neel-ordered domain, i.e., it docs 
not have the quantum mechanical magnitude y S(S + 1) 
which might have been naively expected. The full low-T 
impurity susceptibility is hence xfmp = \t^ 2 + fXimp 


In order to have a simple system which reproduces 
this expected low-T behavior, and also accounts for non- 
universal behavior at higher T for different types of im- 
purities, we introduce a simple effective few-spin model. 
The idea is to model a large Neel-ordered region by a 
classical vector N. Three Hamiltonians corresponding to 
the original models are defined: 

Hf = aS -S e + rN -S e -h z M§, (4a) 
H cB = rN ■ S e - h z M z _, (4b) 
Hf = aS Q ■ S e + rN • S e + J x S a ■ S - h z M%. (4c) 

The effective Hamiltonian for the vacancy model, H_, 
contains a single spin S e representing a remnant spin-^ 
due to the sublattice asymmetry caused by the vacancy. 
It is coupled to the unit vector N, representing the orien- 
tation of the local Neel order of the host antiferromagnet. 
The magnitude of this order is absorbed into the effective 
coupling r. In Hf we "reinsert" the spin So that was 
removed in the vacancy model, and couple it to S e . Fur- 
ther, including the Heisenberg interaction between the 
host spin So and the added spin S a , we arrive at the 
effective Hamiltonian for the added spin model Hf . 

We determine the impurity susceptibilities © of the ef- 
fective models in the same way as for the original models. 
In Eqs. Q we have explicitly indicated how an applied 
external field, which defines the z direction, couples to 
the systems. The magnetization operators are M£ = S%, 
Mg = S z + S z , and M* = S z a + S z + S z . The suscepti- 
bilities for k = 0, — , +, corresponding to Hf , ZP ff , and 
Hf , can be written in the form 

Xfc = ^M = ~x! + ~Xfc, (5) 

where || and _L refer to the directions parallel and per- 
pendicular to the vector N. After straight-forward diag- 
onalization of the effective Hamiltonians in the || basis, 
the two components can be evaluated using 

xl = i((Mj) 2 ), (6a) 

xi = I ' dr(M^(r)Mt(0)). (6b) 
Jo 
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FIG. 1: Impurity susceptibilities for different system sizes L 
with one vacancy. Error bars are smaller than the symbols. 
The inset shows a comparison between L = 64 simulation 
data and the effective model. 

In (J5J an average over all orientations of N relative to the 
fixed z axis has been taken. The coupling of the impurity 
spin to N, therefore, and in accordance with Ref.|5|, leads 
to a classical-like low-T divergence; xt^> "~ * T2 T + §X±> 
where the second term is constant at low T (instead of 
having the r = form )■ 

We do not attempt to derive values for the couplings a 
and r [however, J± — J as in lfTc|l ]. Moreover, at least r 
should in principle have some T dependence. All results 
presented will be for constant a/ J — 2.1 and r/J = 1.75, 
which give a reasonable over-all agreement with the SSE 
calculations. We will show that although the effective 
model is highly simplified, it captures some of the differ- 
ences between the vacancy and the added spin. 

Next, we present the results of the SSE calculations for 
the models Q and compare with the corresponding ef- 
fective models. The impurity susceptibilities of interest, 
Eq. ©, are defined as differences between two extensive 
quantities. Although improved estimators were uti- 
lized in the SSE, the relative statistical noise grows very 
rapidly with increasing L. We are therefore currently 
limited to L < 64 and T > J/32. 

The single-vacancy impurity susceptibility multiplied 
by 4T is shown in Fig. ^ At high T the data for different 
L are indistinguishable, while at low T finite-size effects 
are clearly seen for L < 16. In the range of T considered, 
all finite-size effects are eliminated within statistical er- 
rors for L = 64. The observed behavior at high T is 
easily understood as the total susceptibility is then just 
the sum of Curie contributions of each independent spin; 
xfi p (T - 00) = (L 2 - 1) /AT - L 2 /4T = -if At low 
T the S = I ground state of £f_ and S = of Ho lead 
to the observed Curie behavior Ximp = ^ or smau L. 



3 



1 .0 


v o added spin, L=64 
n£ effective model 


0.8 


gi • vacancy, L=64 
effective model 


0.6 

H 

cN 




i — i 


•a \^ 


7 0.4 




+J 

0.2 




0.0 


— "V - 

«H. • 



0.1 1.0 10.0 

T/J 



FIG. 2: Xi^> ~ T2 7 ^ OI s y s t ems of size L = 64 with a vacancy 
and an added spin. The straight lines are fits to the low-T 
simulation data. 

As L grows this finite-size effect vanishes, and we observe 
the predicted A A contribution [5j arising from the longi- 
tudinal part, in analogy with the effective model. Hence 
the nearly constant 4Txj^p ~ g for L = 32 and 64. The 
effective model reproduces the behavior reasonably well, 
as shown in the inset of Fig. ^ 

Looking more carefully at the SSE data, the predicted 
low-T constant behavior in the quantity x£np ~ TiTi 
which should reduce to fx^J as T — » 0, is not observed. 

fn Fig. (solid circles) we present results for ximp ~ J2T 
for L = 64 (data for smaller L indicate that there are 
no significant finite-size effects for L = 64). ft shows 
an apparent logarithmically divergent behavior, roughly 
from the onset T of renormalized-classical behavior in the 
pure 2D Heisenberg susceptibility, T/J w 0.3 18]. This 
is also approximately where the corresponding effective- 
model result converges to a constant. 

Results for two vacancies are shown in Fig. [3] The 
high-T behavior has the same explanation as for the sin- 
gle vacancy. At low T, the moments due to the two va- 
cancies, which reside on different sublattices, are pinned 
by the Neel order antiparallel to each other, resulting in 
a vanishing xjmp- Hence, in this case Ximp does not di- 
verge for small L and we do not, therefore, multiply our 
results by T. The inset in Fig. shows a comparison 
between one and two vacancies for L = 16. The T at 
which the two curves deviate corresponds to a correla- 
tion length of the same order as the separation between 
the two impurities, £ f=a L/2, i.e., above this T the two 
impurities couple to different Neel domains and behave 
as independent vacancies. Since £ diverges exponentially, 
the point of deviation moves very slowly towards T = 
as L grows. At low T the resulting xtmp ^ or ^ w0 vacanc i es 
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FIG. 3: Impurity susceptibility for systems with two vacan- 
cies. The inset shows results for L — 16 systems with one and 
two vacancies. 



shows little T dependence for large L. However, no sign 
of convergence of the plateau value is seen as the system 
size grows. If xtnp ^ s finite as T — > 0, we would expect 
such a convergence at low T (with a peak at interme- 
diate T for very large L, where the cut-off of the 
divergence occurs at low T) . It should be noted that in a 
finite system there will always be some interactions also 
between the _L components of the two impurity spins at 
low T, We can therefore not expect the behavior for two 
vacancies to be given in a straight-forward way by the 
single- vacancy results in Fig. [21 The roughly In (L) di- 
vergence of the plateau height in Fig. is, however, fully 
in line with a log-divergent Ximp ^ or a sm gle vacancy. 

Results for ximp f° r systems with an off-plane added 
spin are shown in Fig. 0] The high-T behavior, as well 
as the low-T behavior for small systems (L = 4,8), is 
understood by the same arguments as for the vacancy. 
Note that in this case the high-T impurity susceptibil- 
ity is instead of — for the vacancy. Again, 
we believe that all finite-size effects are eliminated for 
L — 64 down to T = J/32. Contrary to the behavior for 
the vacancy in Fig. ^ the T at which 4Txf^, assumes 
an almost constant value | has not yet been reached at 
T = J/32. Moreover, we note the substantial differences 
between yf'"t and yf '~ at intermediate T. In particu- 
lar, the shoulder seen in Fig. 0]for T ~ 0.2 — 1.2 has no 
counterpart in Fig. 2] This feature is clearly due to the 
internal structure of the added spin impurity and is re- 
produced very well by the effective model. In Fig.[2](open 
circles) we present the results for Ximp ~ T2 for T = 64. 
As in the vacancy case, it appears to be log-divergent at 
low T. Note also that the effective model describes the 
added spin impurity very well down to quite low tern- 
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FIG. 4: Impurity susceptibilities for different system sizes L 
with an off-plane added spin. The inset shows a comparison 
with the effective model. 



peratures. It captures the behavior better than for the 
vacancy, but a better agreement for the vacancy can also 
be achieved by using slightly different a and r. 

Summarizing our results, we have confirmed the pre- 
dicted [f| contribution to the impurity susceptibil- 
ity. However, we also find a logarithmic divergence of 
Xlinp — T2T' mstea( l 01 tne constant §X[Jr np predicted for 
this quantity at low T (and reproduced with our effective 
model). We cannot, of course, completely exclude an ap- 
proach to a constant at still lower temperatures, but such 
a slow convergence had also not been anticipated. We 
note that the separation of xf mp into transverse and lon- 
gitudinal components, Eq. (J5J, is strictly correct for the 
2D Hcisenberg model only at T = 0. However, the very 
sudden cutoff of the divergence seen in the two-vacancy 
data in Fig. [21 supports the notion of a component align- 
ing very strongly to the local Neel order (which becomes 
the global order at the L-dependent crossover T) and 

X 



justifies the relation |xfm P 



uup j^^t also at rela- 
tively high temperatures. Hence, our results are most 
naturally interpreted as a log divergent Ximp- This is in 
fact in line with a Green's function calculation by Na- 
gaosa et al. 0. They found that for a system with a 
vacancy, the frequency dependent impurity susceptibil- 
ity at T = 0, Xim P (T = 0, u>), was log divergent when 
u; — s- 0. In view of the renormalized-classical |l2| picture, 
this is consistent with exact results for the classical 
2D Heisenberg model. An anomalous perpendicular sus- 
ceptibility was also recently noted for the 5=1/2 model 
at finite impurity concentration 20]. 

The results presented here call for a reexamination 
of the field theory [f| of quantum impurities in the 
renormalized-classical regime. Very recent efforts to ex- 



plain the log divergence, motivated by our numerical find- 
ings, have indicated that the impurity moment acquires 
a previously unnoticed correction of order Tln(l/T) to 
its lea ding order value S in the renormalized-classical 
regime [2l], |23 , and that this can be interpreted as a log- 
divergent contribution to Xfinp as proposed here (results 
in the quantum-critical regime |f| remain unchanged). 

Apart from the log divergence, the effective few-spin 
model that we have introduced here gives a good de- 
scription of the impurity susceptibility at high and in- 
termediate T. In particular, it captures very well the 
nonmonotonic T dependence that we have found in the 
case of an added spin. 
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